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O "■ Abstract 
O 
O 

We present an entirely analytic, leading log order determination of the fric- 
■ tion an electroweak bubble wall feels during a first order electroweak phase 

transition. The friction is dominated by W bosons, and gives a wall velocity 
parametrically ~ a w , and numerically small, ~ .01 — 0.1 depending on the 
Higgs mass. 

I. INTRODUCTION 
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Electroweak baryogenesis is the name for the production of the baryon number asym- 
metry of the universe at the electroweak epoch. It is possible in extensions of the standard 
^ ' model where the three famous Sakharov conditions can be met: 

^ ' 1. Baryon number violation is efficient, 

q h ' 2. The discrete symmetries C and CP are violated, and 

<D ' 

3. There is a departure from equilibrium, coincident with baryon number violation turn- 
ing off. 
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All three of these criteria potentially exist in the standard model. Baryon number violation 
is not only present, but efficient [|]||. The expansion of the universe provides departure 
from equilibrium. And of course, C and CP are known not to be true symmetries. 

Whether or not electroweak baryogenesis can explain the size of the observed baryon 
number of the universe, which is about || 

baryons _ , n nX _ _ 1 n 



entropy 



2 — 7)xl0" li , (1.1) 



is a more detailed question. In the minimal standard model, both the available CP violation 
|4j], and the departure from equilibrium ||, appear to be grossly insufficient. However the 
question is quite open in extensions of the standard model. 

Baryogenesis at a first order electroweak phase transition is a complicated process, and 
requires understanding several things. First, we must be able to compute the strength of 
the phase transition. The tools for doing this are now well developed |5]|J. Next, we must 
be able to compute the efficiency of baryon number violation; here the tools are also well 
developed |7HT2| . Finally one must be able to compute the microscopic dynamics of baryon 



number carrying excitations in the presence of an electroweak phase interface, henceforth 
called the "bubble wall" |13|] . This problem is not well under control, but the papers quoted 
show that the results strongly depend on another factor, which is the velocity of propagation 
of the bubble wall. This paper will discuss the computation of this bubble wall velocity, 
which is important for baryogenesis, and is also interesting as an example where such a 
dynamical quantity can be computed from first principles. 

A fairly substantial literature already exists on the electroweak bubble wall velocity. 
Since a bubble wall liberates latent heat as it propagates, hydrodynamic considerations are 
potentially important. Hydrodynamic considerations appear in [|Tl|-|rS(]. The conclusion is 
that, if the friction on the bubble wall is small, then the hydrodynamics are important; but 
if the bubble friction is large, so the wall velocity is small, then all that matters is the general 
rise in temperature from the bubbles in aggregate as the transition proceeds. 

There is also a literature on the friction the bubble wall feels |T9|-p5|. The paper of 
Khlebnikov []35] shows how the friction is related to the self-energy of the zero mode of the 
bubble wall. That is, the friction arises from the back-reaction on the wall of the disturbance 
from equilibrium of excitations, induced by the motion of the wall. Most of the papers quoted 
study this back reaction by treating the excitations with kinetic theory. The exception is 
|25| |, where it is argued that infrared SU(2) gauge field and Higgs field excitations are most 
important, and that it is more appropriate to treat them as classical fields, which can be 
done nonperturbatively on the lattice. 

Both the kinetic descriptions, and classical nonperturbative treatment, missed one im- 
portant piece of physics, however, which as we will see leads them to be incorrect paramet- 
rically. That is, they miss the physics of screening and Landau damping, which dominates 
the dynamics of infrared gauge fields. The importance of this physics has been pointed 
out by Arnold, Son, and Yaffe, |]T|, |2l||2"T|] , in the context of determining the baryon number 



violation rate, and has been further discussed in P ^P9[[L2| . The central result is that the 
SU(2) gauge field A, instead of evolving under (classical) equations of motion of the form 
(in temporal gauge, ignoring nonlinearities) 

d ^ = -(k 2 + m 2 )A, (1.2) 

with m 2 = (g 2 /4)cj) 2 the mass squared induced by a Higgs condensate, instead undergoes 
overdamped evolution, 

7 lll3l = -(k 2 + m 2 ) A + noise , (1.3) 
4k dt 

with mj5 the Debye mass squared, = (11/6) g 2 T 2 in the standard model. In this letter 
we will see what consequences this has for the kinetic description of the bubble wall friction. 



II. INGREDIENTS 

For concreteness we will work here with the minimal standard model, even though baryo- 
genesis in that model is ruled out. Since only the gauge fields are overdamped, and since 
only the SU(2) and U(l) fields have large interactions with the Higgs condensate, the ex- 
tension to models with more fields, such as the minimal supersymmetric standard model, 
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should be straightforward. From here on we will also neglect the U(l) field, which is the 
same as taking g' ^ g. This is probably a reasonable approximation. Going beyond it would 
increase the bookkeeping but would not make our treatment substantially more difficult. 

This work will be strictly analytic and strictly based on parametric expansions in formally 
small quantities. In order to obtain an electroweak phase transition of a strength which can 
be analyzed perturbatively, we will take 

3 <A«07log(l/<7), (2.1) 



where A is the Higgs self-coupling in the same normalization as used in pOfl . As we see in a 
moment, the gauge field condensate is parametrically 4> ~ g 3 T/X, so the induced gauge field 
masses are #0o/2 ~ g 4 T/X. The longitudinal gauge fields have a Debye mass ~ gT ^> g 4 T/X 
and can be neglected [pOf| ; the one loop effective potential is then approximately [f2(| 



V, ioo P (0) = - -f^0 3 + \^ . (2.2) 

A broken minimum exists if dV/dcf) = at some nonzero 0o- The broken minimum is 
degenerate with the minimum at <fi = at 

m 2 (T ) = y oi , (2.3) 



and the value of is 



WTJ = g. (2.4) 



The cf) profile of the electroweak bubble wall at T eq is given at leading order by 
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(2.5) 
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Here z is a space coordinate orthogonal to the bubble wall. 
The gauge field mass times the wall width, 




# 16ttV2A g*T ; _, ^ 

Lm w — L— — — ^r^TT = TlT >>1) 2 - 6 
2 g 6 T 167rA o 0o V A 

is large, and it is therefore possible to treat gauge field excitations, in the broken phase and 
inside the bubble wall, with kinetic theory. The kinetic theory description, which amounts 
to taking the Higgs field background as approximately homogeneous and expanding in its 
gradients, breaks down at the small z (symmetric phase) edge of the bubble wall, when 

(<p/(po) ~ \J^J 9 2 - 111 fact, since A = (g 2 /S) is the condition for ran = raw-, the kinetic 
description works fairly well at remarkably large Higgs masses. 

All of the equations above are valid at leading order in X/g 2 or g 3 /X. when the former 
breaks down, higher loop corrections and corrections from Higgs loops (neglected here) 
become important. When the latter breaks down, the longitudinal gauge fields become 
important. We also need g 3 A to ensure that the relevant gauge fields will be overdamped. 
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III. FRICTION 



The friction on an electroweak bubble wall is defined as the excess pressure on the wall 
(directed towards the broken phase), over the equilibrium value; 

friction = P - P cq = P + V((f) = 0) - V((f> ) = P + AV . (3.1) 

We expect the friction to depend on the bubble wall velocity, with the condition P = 
determining the steady state bubble wall velocity v w , which is what we want to know. We 
can define a linear response friction coefficient rj as a limit 

,. AV steady state friction 
77= hm = . 3.2 

AV^O Vw V w 

This is what we want to determine. 

The friction on the bubble wall depends on the departure from equilibrium of the plasma 
excitations inside the bubble wall. In a kinetic theory description, the excitations are de- 
scribed by population functions f(k,x). We write them as an equilibrium part /o, 



f ° = e M E/T)±l ' B =l/*" + ">'<*>. < 3 ' 3 > 

with + for fermions and — for bosons, which is the case we will care about, plus a departure 
from equilibrium Sf. The friction a bubble wall feels is, in the kinetic description |l~9|- [2T|j24]] p|, 



friction 



DOF 



f°° r f d?k d 4> dm 2 dE . s 

L dz LJwrTz^d^ 6f ' (3 - 4) 



where the sum is over degrees of freedom which get a mass from the Higgs field. This 
expression has a clear intuitive meaning; it is the sum over excess particles of (dE/dz), the 
force the wall exerts on them. It remains to determine Sf and evaluate the integral. 

The friction will be dominated by the gauge boson contribution. The contribution from 
Higgs bosons is smaller because their dm 2 /dip is smaller by X/g 2 , and because their evolution 
is not overdamped. The fermionic contributions, such as that from the top quark, are smaller 
because Fermi-Dirac statistics lack the infrared divergence of Bose-Einstein statistics, so 
their 5f has much weaker infrared behavior. The parametric argument appears in ||25|| , 



and a reasonable estimate of friction from top quarks appears in [2l|; it proves numerically 



smaller than what we find below. The case where there is a light scalar top is more difficult 
and we do not consider it; for the MSSM the friction we find should be viewed as a lower 
bound rather than a tight estimate. 

As discussed above, the gauge fields undergo overdamped evolution given by Eq. fll.3|). 
Since / oc A 2 , the equation for / is 

k™>d df 2 . 

— - = -Ef + noue, (3.5) 



L A derivation of sorts can be found in [E4|, but the expression is implicit in the earlier references as well. 
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where the noise is of the right size to ensure that, for m 2 time independent, / will approach 
f ; so averaging over the noise, 

(k 2 + m 2 )f + noise -> E 2 5f . (3.6) 

Also, df /dt = dfo/dt + d(8f)/dt. At small v w , the limit we are interested in, 5f <C f , and 
d(5f) I dt may be dropped. Further, 

dfo _ dcfi dm 2 dE df _ d<p dm 2 1 . 

dt ~ dt d0 dm 2 d£ ~ ^ d0 2 £T M +/oj ' 1 j 

Therefore, the departure from equilibrium of a gauge boson degree of freedom is 

7rm|^ d<p dm 2 w 

Note that transport plays no role in setting 5f; this is because the gauge fields are over- 
damped. Substituting this into Eq. (3~3), and noting that there are 6 species of transverse 
W bosons (3 flavors times 2 spins), gives 

QTTv w ml r°° (d(f)dm 2 w \ 2 f d 3 k /„(l + / ) 
faction = dz[ — —— / . o — ^- . (3.9) 



-oo \dz d(j) J J (2tt) 3 AkE A T 

Since /o is monotonically decreasing, the momentum integral is infrared dominated, cut off 
by the nonvanishing W boson mass. Therefore it is appropriate to make the approximation, 
valid in the infrared, that /q ~ 1 + fo ~ T/E, and evaluate the integral; 



d 3 k / (l + /o) /•°°A; 2 dA; T T /■«> d(E 2 ) T 



(2tt) 3 4fc£ 4 T Jo 2vr 2 4kE e 16tt 2 ^ -E 6 327r 2 m^ 
The friction, using mw = g<p/2, is then 



(3.10) 



w 



3mlT r°° , /d0\ 2 1 



friction = v w —- — / dz [ — ) — . (3.11) 

327r J-oo \dz (b z 



(3.12) 



Writing (d<j)/dz)dz = d<p, and using Eq. (|2.5|) to write 

d0 20(0 O - 0) 



cfe L0 O 



Eq. ( |3.11|) now gives 



= Sm^T x / /*> (0Q - 0)d0 = /* (1 - x)dx \ ^ 



167rL yJo 0o0 •/() x 

There is a log divergence arising from the symmetric phase side of the bubble wall. The log 
will be cut off where the first approximation used to derive Eq. ( |3.11| ) breaks down. The 
perturbative expansion breaks down when my/ ~ 9 2 T, or ~ gT, which is at (0/0o) = 
(\/g 2 ). The kinetic theory description breaks down at (0/0o) = \J~\J~g 2 ~, see Eq. (|2.6| ); this 
occurs first. 
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Since the degrees of freedom which dominate the friction are those with k ~ m, when m 
drops below 1/L, it is no longer appropriate to treat the particles as seeing a slowly varying 
wall. Such degrees of freedom see a wall which is sharper than their wavelength can resolve. 
For those degrees of freedom with kL ^> 1, we find the friction scales as 1/L. This must go 
over to an L independent value for wavelengths which cannot resolve the thickness of the 
wall, which means that their contribution is less than the kinetic theory estimate. Hence 
the log is cut off at ((j)/4>o) • 



'\/g 2 , the contribution from very infrared degrees of freedom 
is sub dominant. Hence the friction we determine is 



77 



3mlT 
16nL 



\og{m w L)+0{l) 



3 m D 9 2t 4 [, 9 , nf.s 



(3.14) 



The first expression makes no assumptions about the effective potential or wall thickness 
and should be valid in extensions as well as the standard model. 

Now we comment on the parametric form of rj. Taking A ~ g 2 and neglecting logs, the 

We can guess this on dimensional grounds by noting that 



friction coefficient is rj oc a^T 



pressure 




energy x time 


velocity 




length 4 



(3.15) 



and that rj arises from infrared gauge field physics. Such physics has a natural energy 
scale ~ T, a natural length scale ~ l/a w T, and a natural time scale ~ l/(a^T) f7]; so on 
dimensional grounds we should have anticipated 77 ~ &wT A . This is to be contrasted with 
the pressure driving the bubble wall, which by the same parametric estimates must be 

P~a 3 T 4 . (3.16) 



P 



energy 



length 3 

Hence, the bubble wall velocity is parametrically v w ~ a 



IV. BODEKER'S EFFECTIVE THEORY 

The reason we considered the case A g 2 / log(l/g), rather than A g 2 , is because in 
the parametric regime 



a 2 

< A < g 2 (4.1) 



log(l/») 

The gauge bosons with m ~ g0o/2 do not obey Eq. ( |1.3j ); instead Bodeker's effective theory 
is applicable 0; 

,_ = - B M + no 1S e, g = 3g2r|og ° /g) . (4.2) 

up to corrections suppressed by log(l/g). In this case the derivation proceeds analogously, 
but the behavior is slightly less infrared dominated, and no log occurs in the integral. The 
final expression is 



1 



256^ log(l/#) \y/Xj \g 2 T 



2 / w 



a 2 w T 4 . (4.3) 
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This result is in parametric agreement with what we found above. It is only applicable in 
an extremely narrow parametric range, and receives corrections suppressed only by {X/g 2 ) 
(from the loop expansion) and g 2 log(l/g)/A (from the breakdown of Bodeker's effective 
theory) . 



V. CONCLUSION: VALUE OF THE WALL VELOCITY 

Since we only have the friction at leading log, which means with at least a factor of 2 error, 
we will do with a fairly crude estimate of the pressure which drives the bubble. In references 
20| , |2T|| the nucleation temperature is estimated as occurring at m 2 (T nuc ) = 0.8m 2 (T eq ). For 



this value, using Eqs. ( |2.2|) and (|2.3|) , the pressure driving the bubble wall is approximately 
P = .94 x (# 6 /i(D247rA 3 )c^T 4 , anc i the wall velocity is about 

\ 5 

9 \ a w 




which, for m H ~ m w or A = g 2 /8, and estimating \og(g/^/X) ~ 1, is v w ~ a w /A < 0.01. For 
X/g 2 ~ .04, as required to give a sufficiently strong phase transition, v w ~ 0.1. The numerical 
value for the bubble wall velocity is expected to be small. The numerical estimate from the 
result determined using Bodeker's effective theory is similar; v w ~ 0.2 for X/g 2 = .04. Neither 
estimate is very reliable because the constant under the log has not been determined; but 
the conclusion v w <C 1 is clear. 

The friction we find is larger, and the value of v w smaller, than in previous literature. In 



particular we find a much larger friction than the numerical results of p5[ indicate. This is 
partly because there we did not include hard thermal loop effects for the gauge fields, but 
overly aggressive data fitting may have contributed. It would be interesting to make a new 



numerical analysis using the techniques developed in p"T|. 

We should also briefly comment on how the result may change in extensions to the 
standard model and beyond leading order in v w . Beyond first order in v w , a few effects 
become important. We may not be able to neglect 5f next to f in getting Eq. (|3.8| ); 5f, 
and the friction, will be larger. Also, the frictive pressure will change the bubble wall shape. 
Since most of the friction is on the symmetric phase side of the wall, the wall will become 
narrower, see p4j| , which also increases the friction. Finally, in extensions to the standard 
model (which are the only viable candidates for baryogenesis because they can provide both 
a strong phase transition and the heavy Higgs boson required by experiment), the bubble 
wall is typically thinner, because the Higgs mass is larger. This increases the friction, even 
before considering new contributions from extra light bosons such as a light scalar top. For 
these reasons we anticipate that electroweak bubble wall velocities are quite generally much 
less than 1. 
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